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*^ I A family of explicit exact solutions of Einstein's equations in four and higher dimensions is 

studied which describes the gravitational field of an object accelerating due to an anisotropic 
^ . emission of photons. It is possible to prescribe an arbitrary motion, so that the acceleration 

qh| of such photon rocket need not be uniform — both its magnitude and direction may vary 

with time. Except at location of the rocket the spacetimes have no curvature singularities, 
and topological defects like cosmic strings are also absent. Any value of a cosmological 
constant is allowed. We investigate some particular examples of motion, namely a straight 
CN ■ flight and a circular trajectory, and we derive the corresponding radiation patterns and the 

mass loss of the rockets. We also demonstrate the absence of "gravitational aberration" in 
00 '. such spacetimes. This interesting member of the higher-dimensional Robinson-Trautman 

' class of pure radiation spacetimes of algebraic type D generalises the class of Kinnersley's 

solutions that has long been known in four-dimensional general relativity. 
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1 Introduction 

X 

^ ' In 1969, Kinnersley [1] introduced and investigated a class of exact spacetimes that describe 

a localised object which accelerates due to the back reaction of the emitted null radiation. 
In the axially symmetric case the object moves along a straight line corresponding to the 
axis, but more general solutions also exist in which the motion is completely arbitrary. They 
may serve as interesting self-consistent exact models for the motion a photon rocket that 
is propelled by a specific anisotropic emission of photons. Such solutions have attracted 
considerable attention pl-fTO] , in particular due to some unusual properties of the associated 
gravitational radiation. 

The Kinnersley solutions are of algebraic type D and belong to a large family of Robinson- 
Trautman spacetimes [llH14j which is defined by the property that it admits a geodesic, 
shear-free, twist-free but expanding null vector field. Interestingly, within the class of 
Robinson-Trautman spacetimes there are other distinct type D solutions which also de- 
scribe accelerated sources, namely the famous C-metric. This is a specific vacuum solution 
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which represents a pair of black holes, uniformly accelerating under the influence of cosmic 
strings or struts along the axis of symmetry p3] (for a recent review see [Hj). In contrast 
to the C-metric, the spacetimes describing Kinnersley's photon rockets are non-vacuum (the 
rockets emit pure radiation) and, apart from the location of the rocket, they are regular 
everywhere. 

It is desirable to find and study higher-dimensional extensions of such classes of exact 
spacetimes which represent accelerated objects. There have been many attempts to generalise 
the C-metric to higher dimensions, but it has not (yet) been found. Interestingly enough, an 
exact class of pure radiation spacetimes which includes those of the Kinnersley rockets in an 
arbitrary dimension has been recently discovered within the Kerr-Schild family |16yi7) and 
independently in the Robinson-Trautman family }18]. 

It is the purpose of this contribution to present and analyse such spacetimes. In sec- 
tion [2] we briefly summarise the Kinnersley photon rockets in four dimensions. Subsequently, 
in section [3] we introduce very useful Newman-Unti coordinates adapted to an arbitrarily 
moving test particle in a flat space of dimension D. In section H] we present a generali- 
sation of the Kinnersley solution to higher dimensions, discussing the Kerr-Schild and the 
Robinson-Trautman forms and a possible cosmological constant. Sections [5] and [6] concen- 
trate on particular situations in which the rockets accelerate along a single spatial direction 
or move along a circular trajectory. In both these cases we derive and plot the corresponding 
radiation pattern and the mass loss formula. In the final section [7] we give the Christoffel 
symbols necessary for study of geodesies in these spacetimes and we demonstrate the absence 
of "gravitational aberration". 



2 Kinnersley photon rockets in = 4 

Kinnersley's solution [IJ belongs to the family of Robinson-Trautman space-times with an 
aligned pure radiation. In four dimensions, the metric can thus be written in the standard 
form dMl] 

ds^ = 2^dQdC-2dudr -2Hdv?, (1) 

where 

2H = K-2r{\ogP),^-'^-^-K\ (2) 

V o 

Here A is the cosmological constant while P is an arbitrary function P{u,C,C)- The func- 
tion K{u, C, C) = A log P, where A = 2P^d(^d^ , determines the Gaussian curvature of the 
2-surfaces spanned by the (complex) spatial coordinate on which r = 1 and u is any con- 
stant. 

The family of Kinnersley rockets is obtained when 

P = A{u) + B{u)(: + B{u)C + C{u)(:C, (3) 

in which A, B, C are functions of u (B may be complex), so that 

K{u) = 2{AC - BB) . (4) 

For > 0, these solutions represent the gravitational field of an arbitrarily moving object 
located at r = whose velocity is encoded in the functions A,B,C (for the explicit relations 
see equation ([28|) below; more details are given in section [3]). Since the Gaussian curvature 
K{u) is independent of the transverse spatial coordinates there are no poles over the 
compact surfaces u = const, (at any r) which means that for r ^ these spacetimes are 
everywhere regular, free of curvature singularities and cosmic strings or struts. Indeed, 
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using the natural null tetrad k = dr, I = du — Hdr, m = (P/r) the only component of the 
curvature tensor for the metric ([T]) with ([3D is 

m(u) 

^2 = -^. (5) 

The spacetimes are thus of algebraic type D, and there is a curvature singularity only at 
r = where ^2 diverges. Conformal infinity is located at r = 00 where the space-times 
asymptotically become conformally flat and also vacuum, i.e., Minkowski or (anti-)de Sitter, 
according to the sign of A. 

The spacetimes contain an aligned pure radiation field (that is flow of matter of zero 
rest-mass, emitted from the source located at r = 0) with an energy-momentum tensor of 
the form T^;^ = pk^k^ . The radiation density is 

P = -2 ' (6) 

where the function n? is determined by Einstein's equation as 

47rn^ = 3m(logP),t, . (7) 

In the particular case when 5 = and K = +1, the Kinnersley photon rocket moves 
along a single axis. Performing a suitable transformation, the metric ([I])-(l3|) then can be put 
into the form 

ds =— 1 r — 2arcosw — a r sin v \a.u 

- 2 dn dr 2a sin ^<lud^+ (di?^ + sin^ i9 d(/)^) , (8) 

in which a{u) is an arbitrary "acceleration function" of the coordinate u, see [T][2|llUj for more 
details. This represents a singular source located at r = 0, of decreasing mass determined by 
m{u), which emits pure radiation and accelerates in Minkowski, de Sitter or anti-de Sitter 
space along a straight line (which is the axis of symmetry) due to the corresponding net back 
reaction. It thus serves as a simple exact model of a rocket that is propelled by the anisotropic 
emission of photons whose radiation field profile is determined by the field equation ([7]) as 
n^(n, -d) = -^[—m^uiu) + 3 a{u) m{u) cos -d]. 



3 Arbitrarily accelerated coordinates in flat space 

Before presenting a generalisation of the Kinnersley photon rockets solution to higher di- 
mensions (in the following section 2]), it will be important to introduce and summarise a 
particularly useful coordinate system adapted to an arbitrarily moving test particle in a flat 
space of dimension D. For D = A this was described already by Newman and Unti in their 
classic work [19j, and subsequently elsewhere. 

Let us consider Minkowski space of dimension D with the standard metric 

ds^ = r/„/3 dZ"dZ^ = -(dZO)^ + {dZ^f + ■■■ + (dZ^^l)^ (9) 
and any timelike worldline 

z%u). (10) 

This describes motion of a test particle in the flat space @ , and the parameter u is assumed 
to be its proper time. The corresponding velocity is 

dz" d 
u = z'^da, where i"(ti) = — and da = ■ (H) 
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At any event V = {z^{u), . . . , z^^^{u)} on the trajectory (jlOp we now take the future null 
cone and we label this hypersurface as n = const, which is exactly the value of the proper 
time of the particle when it passes through the vertex point V, as shown in figure [TJ 




Figure 1: Construction of the coordinate system that is naturally adapted to an arbitrarily moving 
test particle in D-dimensional flat space (spatial coordinates are suppressed). At 

any event V of the particle's worldline with velocity w, the future null cone is labeled as u and it 
is assigned the value of the corresponding proper time. The coordinate r is an affine parameter 
along null geodesies which connect V with any event Q, on the cone. These are generated by 
properly normalised null vectors h. 

Each of such null cones is generated by a family of null vectors k = da along the particle's 
trajectory. It is convenient to prescribe the normalisation condition k ■ u = —1, in addition 
to those already mentioned, u ■ u = —1 and k ■ k = 0, i.e.: 




The background coordinates of any event Q = {Z^, . . . , Z^~^} on the null cone u = const, 
can obviously be expressed as = z°'{u) + rk°', where k" represents the generator of the 
null geodesies connecting Q with V, and r is a new coordinate. In fact, r is the correspond- 
ing affine parameter since dr = Z"^ da = k'^da = k, and using the relation (jl2p this can be 
explicitly expressed as r = —ria[s{Z'^ — z"(u)) z^ . Now it only remains to parameterise the 
null vectors k. In the natural frame of the background coordinates of Q, without loss of 
generality all these vectors can be conveniently represented as 




D-l 



(15) 



i.e., k^ = p ^,k^ =p ^n*. Here p{u,n) is some function, and the parameters n' can be 
understood as components of a spatial unit vector n = n^di. Indeed, 

D-l 
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there are D — 1 spatial parameters n*, but since they are constrained by the normaUsation 
condition (|16|) . only D — 2 of them are independent. These new coordinates are adapted 
to an arbitrarily moving (accelerating) test particle which is always located at r = 0, the 
vertex of all future-oriented null cones. They are labelled by the coordinate u which has the 
meaning of a retarded time because it coincides with the proper time of the moving particle. 

Now it is straightforward to express the Minkowski space Q in terms of these new 
coordinates {u,r, n*} using the relation 

Z" = z°(n)+rA;"(u,n). (17) 

Employing relations (fT2]) . (fT3]l . (fH|) and also the identities 

Irj^^k^dk^ = dirjapk^'k^) = 0, (18) 
r]^pi''dk'^ = (logp),,du, (19) 

r]ai3 dk^'dk'^ = —Y. d"*^^'' ' (20) 



we arrive at the following metric form 

„2 D-l 



ds^ = ^ ^ 5ij dn'dn^ - 2dudr - (l - 2r(logp),„) du^ , (21) 



?.,j=i 

with the constraint ()16p . The function p(u,n) in ()2ip is not arbitrary since the condition 
k ■ u = —1, i.e. p^ . must be satisfied. It immediately follows from and (fT5|) that 

D-l 

p(n, n) = i°(n) - ^ n'z\u) , (22) 

where z"{u) are components of the velocity (|lip of an arbitrarily moving particle whose 
worldline is z°'{u). 

3.1 Flat metric in accelerated spherical-like coordinates 

The most natural explicit parametrisation of the components n* of the unit vector n, which 
satisfies (|16p . is given hy D — 2 spherical angles. In = 4 the standard choice 

= cos 6 , r? = sin Q cos , r? = sin Q sin (/> , (23) 

brings the metric ()2ip to the form 

2 

ds2 = I- (d6»2 + sin^ ^d02)-2dndr- (l - 2r( logp),„) dn^ (24) 

where, due to ([22]) . 

p(ii, 0, (/)) = i*^ — i""^ cos 9 — sin cos </> — sin sin </> . (25) 
This can further be rewritten in terms of the complex stereographic coordinate C ■, 

C = \/2 tan ^ e'^ , (26) 

as 



ds^ = 2 dC dC - 2 du dr - [1 - 2r{ log P),„ j d-u^ (27) 



where the function P = (1 + hCOP is given by 



P(n,C,C) = (i° - i^) - ^(i' - ii')C - ^(i' + ii3)C + + ^')CC. (28) 

This Robinson-Trautman form ([1]), (l2|) of the flat metric, which was presented e.g. in [81119^ 
|20j, gives an exphcit physical meaning to the mathematical coefficients A(u) , B (u) , C {u) in 
the function Notice also that K = 1 due to ([IB]) , as required by equation 

Spherical parametrisation of the components of a unit vector n by D — 2 angles 6i can 
be given in any dimension D. Let us consider 

i-l 

= cos6'i [|sin6'j , (29) 
i=i 

and define 



00-1=0, Y[shiej = l. (30) 

i=i 

Relation (fT6]l is identically satisfied and the spatial part of the metric ([2T]) becomes 

D-l D-2 i-l 



Y,Sijdri'dn^ = ^ (]^sin2%)d6',2 (31) 

i=i j=i 

del + sin^ 01 d^i + sin^ Oi sin^ ^2 d^f + • • • + sin^ 0^ . . . sin^ ^15-3 d0|)_2 , 



where 6*1, 6*2, • • • , 6*0-3 £ [0, vr], 0d-2 G [0, 27r). Clearly, it has the geometry of a sphere ^. 
Flat space in arbitrarily accelerated coordinates thus takes the metric form 



^2 D-2 ^ i-l 



ds^ = ^ ^ ( n sin^ Oj'^def - 2d'udr - (l - 2r(logp),„) du^, (32) 

in which 



P i=i j=i 



D-l i-l 

p{u, Oi) = i°(M) - J2 ^*(^) cos 0i n sin % . (33) 
i=i j=i 

Recall again that the functions z^{u), z^{u), . . . , z^^^{u) are components of the velocity pT]) 
of the test particle located at r = which moves along an arbitrary worldline (jlOp . and u is 
its proper time. 

3.2 Flat metric in accelerated Cartesian-like coordinates 



Another important parametrisation of the general metric (j2ip is obtained in terms of Cartesian- 
type coordinates. These are introduced by 



= TT' i,j,k,l = l,2,...,D-2, 



^i 

1 6ij X X"^ 

— A7 I T*^ T*^ 



^D-1 ^ (34) 
1 + ihi x^x' 



explicitly conformally flat. Indeed, J2i'j=i Sijn^n^ = 1 and 



(with the summation convention over i,j, k, I) in which the spatial part of the metric becomes 

D-l 

(l + i4zxV)2 
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The metric (1211) thus becomes 



2 

ds^ = 5ij dx'dx^ - 2 dii dr - (l - 2r( log du^ (36) 
in which the function 

P = + (37) 

with p given by ()22p takes the form 

P = (iO - i^-i) - (5i^. ii) + i(iO + i^-i) x^x^' . (38) 

This is exactly the canonical form of the D-dimensional Robinson-Trautman metric, as given 
by equations (37), (38) of [18] in the flat-space case when A = and fi = 0. In fact, it gives 
an explicit physical meaning to general coefficients in the function (42) therein, 

P = A{u) + Bi{u) x' + C{u) b^j x'x^ , (39) 

namely 

A{u) = z^{u)-i^-\u), 

B,{u) = -i\u), (40) 

w 



C{u) EE l[z\u) + z''-\u)]. 



Consequently, the Ricci scalar IZ corresponding to the spatial metric hij = P ^6ij is 



D-2 



^^""^ -AAC-J2Bf = -Vapz-z^ = l. (41) 



{D-2)iD-3) 



i=l 



The transverse (L* — 2)-dimensional space, covered by the coordinates x*, thus has a constant 
positive curvature, i.e., it is a sphere 5^"^ of radius r, as in p2]) . 

In particular, when the test particle located at the origin r = of the coordinates of (|36|) is 
at rest with respect to the Minkowski background ([9]) then z^{u) = 1, z^{u) = = z^~^{u). 
The metric function (j39|) thus reduces to P = 1 + j 6ij x'^x^ , which gives exactly the well- 
known form of the Schwarzschild-Tangherlini black hole spacetime in the limit when its 
mass parameter vanishes. 

The function ()28p obviously equals ()38p for = 4 if we introduce the complex coordinate 
C, = + i^;^) and relabel z-*^ — t- z'^ — t- — t- z^ 



4 Photon rockets in higher dimensions 

We will now present the class of non-flat solutions which generalise the Kinnersley spacetimes 
from D = 4 to D > 4:. These contain null radiation (representing the emission of photons) 
and a possible cosmological constant A. Such solutions describe an arbitrarily moving photon 
rocket in any dimension. To this end, the coordinate systems naturally adapted to general 
worldlines of test particles in flat space, as described in previous section [3l will be employed 
for the background. 

4.1 The Kerr— Schild form 

In fact, these exact solutions of Einstein's field equations in higher dimensions have already 
been found by several authors. In particular, Giirses and Sarioglu in [16l[T7] presented an 
accelerated metric which has the Kerr-Schild form 

ds^ = 77„^ dZ°dZ^ + 2Vkakp dZ^dZ^^, (42) 
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where the first part represents just Minkowski flat space ([9]), A; is the null vector defined in 
section [31 and the function 2V is 

2m(u) 2A n 

The function m{u) corresponds to the decreasing mass of the accelerating object located at 
r = 0. The associated pure radiation field has the form 

("U, Tl) 

T^u = pk^k^ where p = — , (44) 



and 



„2 ^ + (D - l)m(logP),«. (45) 



D-2 

This is obviously a generalisation of expressions ([7]) valid for = 4. 



4.2 The Robinson— Traut man form 

The same class of exact solutions was independently rediscovered in the context of higher- 
dimensional Robinson-Trautman spacetimes by Podolsky and Ortaggio [18j. It was demon- 
strated that in this class the only solutions with aligned pure radiation and A in any dimension 
D > A are of algebraic type type D and have the form 

ds^ = 7ij dxMx^ - 2 dn dr - dii^ (46) 

where the function 2H is 

T^iu) 2m(u) 2A ^ 

The unimodular metric 7ij(a;^) and the function P{u,x^) must satisfy the field equations 
for the transverse spatial metric hij = P~'^jij, namely that at any section u = uq = const, 
the metric hij{uo,x'') describes an Einstein space, TZij = jj^^TZhij, where TZij and Tl is the 
corresponding spatial Ricci tensor and scalar, respectively. 

The simplest family of such spacetimes arises when the metric hij is of constant curvature 
and thus conformally flat (this is always true when D = 5 since hij is then three-dimensional) . 
In such a case, in suitable coordinates x*, 

7ij = 5ij , (48) 

and, by integrating the remaining field equations, the function P{u, x') must have the general 
form dM]), P = A{u) + Bi{u) x* + C{u) 5ij x'^xK In the fiat-space limit (m 0, A 0) the 
functions A{u), Bi(u), C{u) are related to the components of the velocity it of a test particle 
at r = 0, as given explicitly by expressions ()iO]) . Because this physical interpretation follows 
from relation ()22p . which is the consequence of ()12p . ()15p . and since 



,p k'^z^ = {rja^ + 2Vk^kp) = r,^p k^i^ = -1, (49) 



such an interpretation of the functions A{u), Bi{u),C{u) remains valid also in the general 
(non-flat) case m ^ 0. Moreover, due to the coordinate freedom of the metric ()46p . ()47p given 
hy u = u{u), r = f/u'{u), where u' = g^, which implies P = Pu', TZ = TZu'"^, m = mu'^~^, 
fi^ = -n? u'^ and thus A = Au', Bi = Biu', C = Cu', the corresponding (positive) Ricci 
scalar term in (j47p can always be set equal to 1, see ()4ip . Therefore, the transverse space 



8 



spanned by is a sphere S^^"^ of constant positive curvature. The complete form of such 
Robinson-Trautman pure radiation spacetimes (|46| ) - (j48|) . with m and A non-trivial, thus 
reads 

r2 

ds^ = ^2 ^ij dxMx-^ — 2 dn dr 

-('-^-'■°^n.-fS^- („_,;^„_.) ^-)d.', (50) 

where the function P is explicitly given by ()38p . 

P = (iO - i^-i) - ii) + + i^-i) . (51) 

Here z^{u), z^{u), . . . , i^^^(u) are components of the velocity (jlip of an arbitrarily moving 
photon rocket located at r = 0, whose worldline in the background space ([9]) is z^^u). In 
particular, for an object at rest, P reduces to P = 1 + | Sij x^x^ and metric ([50|) is the D > A 
counterpart of the well-known Vaidya-(anti-)de Sitter metric |131I14] . Such spacetime is 
spherically symmetric and radiation of photons is isotropic with n^(n) = — s— TTi w For 
constant m we recover vacuum Schwarzschild-Kottler-Tangherlini black hole solution. 

It may easily be observed that the metric (|50p can naturally be decomposed into the flat- 
space metric (|36p in the Cartesian-like accelerated coordinates, and the Kerr-Schild term 
2Vdu'^, where the function 2V is given by (flSj) . Indeed, for the null vector k = k°'da = dr 
in the metric (jSOp we obtain ka dZ" = —du. This explicitly demonstrates the complete 
equivalence of the Kerr-Schild form ()42p . ()43p and the Robinson-Trautman form ()50p . ()5ip 
of the metric a photon rocket moving arbitrarily in any dimension D > 4. 

Notice also that the replacement of the flat background metric tj^^ by the curved exact 
metric (7q,/3 = Va/B + 'iVkakp of the photon-rocket spacetime [ ([12]) or (jSUP ] preserves the nor- 
malisations (|12p and (|14p . namely k ■ u = g^p k'^z^ = —1 and k ■ k = g^/B k"k^ = 0. On the 
other hand, the relation (jlSp changes io u ■ u = g^p z'^z^ = — (1 — 2V), where the function 
2V{r,u) is given by expression (|43p . The parameter u thus loses its direct physical mean- 
ing as the proper time of an accelerating test particle located at r = 0. Of course, this is 
not surprising since r = in the complete metric (7^/3 corresponds to a curvature singular- 
ity, and the redshift factor Vl — 21/ = ^1 — 2mr^~^ — (^jj_2^j~)^ r^, which describes the 
time dilation within a gravitational well of the massive photon rocket, also has to be taken 
into account. In fact, similarly as in figure [H the curved spacetime with the metric gap is 
foliated by a family of null hypersurfaces u = const., but the coordinate u now plays the 
role of the time measured by distant observers in asymptotically flat regions r — )• oo where 
y — )• (assuming A = 0; for the case A 7^ see the following subsection 14. 3p . This justifies 
the physical interpretation of the functions i"(n) in ()5ip as components of the velocity u 
of an accelerating rocket, measured with respect to the background flat space (associated 
with the asymptotic regions far away from the rocket). Moreover, the same interpretation is 
obtained in a weak-field limit when the mass of the rocket becomes negligible (m — )• 0). 

Finally, it is possible to employ an alternative spherical-like representation of the ac- 
celerated background coordinates, namely the flat-space metric form ()32p . ()33p given in 
subsection 13.11 The complete spacetime metric gap describing a rocket then becomes 

2 D-2 i-l 

ds^ = ^ 5] ( n sin^ ^j) • - 2 dn dr 
P i=i j=i 

- - 2Klogp),„ - ^ - _ ,)(^ _ 1) r-) dn^ (52) 

where 

D-l i-l 

p(u, 9i) = i°(n) - J2 ^'(^) cos Oi Y[ sin , (53) 
i=i j=i 
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which is fuhy equivalent to the Cartesian-hke form ()50p . ()5ip . 



4.3 The inclusion of a cosmological constant 

In the above metrics, the cosmological constant A may take an arbitrary value. Therefore, 
the solutions represent a photon rocket moving not only in a D-dimensional (asymptoti- 
cally) Minkowski space (when A = 0), but also in de Sitter (A > 0) or anti-de Sitter universe 
(A < 0). These are the three maximally symmetric, conformally flat spacetimes of constant 
curvature R = -J^A. 

For A 7^ 0, the functions z^{u), z^{u), . . . , z^~^{u) in (|5ip or (|5c{p still retain their phys- 
ical meaning as components of the velocity of the photon rocket with respect to the "back- 
ground" frame da corresponding to the coordinates of Q- However, these now need to 
be understood as coordinates of an "external" flat space of dimension D + 1, into which the 
(anti-)de Sitter universe is embedded. 

To be specific, it is well known that the de Sitter and anti-de Sitter spacetimes can be 
represented as a Z)-dimensional hyperboloid 

- {Z'f + {Z'f + ... + {Z^-'f + e {Z^f = ~ '^If ~ , (54) 
embedded in a (D + l)-dimensional flat space with the metric 

ds^ = -{dZ^f + {dZ^f + ■■■ + {dZ'^-^f + e {dZ^f, (55) 

where e = sign A . 

As in (|lUp . we may now consider a timelike worldline z"(n) in this {D + l)-dimensional 
flat space. If we assume that the functions z"{u) satisfy the constraint 

- [z^{n)f + [z\u)f + . . . + [z^-Hu)f + e [z^iu)f = ' " (56) 

at any time u, the photon rocket (with a negligible mass m) during its motion will always 
remain on the hyperboloid ()54p . i.e., in the D-dimensional (anti-)de Sitter universe. The 
corresponding (D + l)-velocity u = i°(u) 9„ is normalised to 

- [z'{u)f + [i\u)]' + ■■■ + [i^-i(n)]2 + e [z'^iu)]' = -1 , (57) 

(for A < we simply modify the metric 7]^^ to diag(— 1, • • • , +1, —1)). This constraint 
may be considered as the relation which determines the value of the function z^{u) in terms 
of the velocity components z^{u), z^{u), . . . , z^^^{u) which occur in the functions ([5T]) or 
(|53p . Effectively, it thus only remains to choose the trajectory of the rocket in such a way 
that it satisfles the condition ([56]) . 

Notice flnally that if the mass of the rocket vanishes, the metric (|5Up or (|52p with m = 
and A 7^ is just the de Sitter or anti-de Sitter D-dimensional space, expressed in coordinates 
whose origin r = is arbitrarily accelerating. An explicit representation of the corresponding 
hyperboloids ()54p in terms of these coordinates can be obtained using the transformations 
(fTTI) and (Uni) as 

70 _ ^0 



Z^ = z"(u) + 



p{u,n) ' 



Z^ = z\u) + - -n\ i = l,...,D, (58) 

p{u,n) 

in which either p4p or ()29p is employed to parameterise the components n* of the unit vector 
n. For the case D = 4 such relations were recently discussed in detail in |10]. 
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5 Motion of the rocket in a single direction 

Let us now concentrate on a particular case in which the photon rocket in D dimensions 
accelerates always in the same spatial direction, i.e., without performing any manoeuvres. 
It may only change its speed. In such a situation, without loss of generality it is possible to 
rotate the spatial axes of the background flat space so that the timelike worldline (jlOp of the 
rocket is 

= {z°(n),z^(u),0,0,...,0}. (59) 
The corresponding velocity ([IT]) \s u = {z^ , i^, 0, 0, . . . , 0} and the condition (fTBj) reduces to 

(i0)2-(il)2 = l, (60) 

The metric of the photon rocket, expressed in the spherical-like coordinates ()52p . ()53p . thus 
becomes 

^2 D-2 ^ i-l 



P' i=i j=i 



2m(u) 2A 9\ , 9 . s 

1 - 2r(logp),„ - ^ - r') dn^ (61) 



where 



p(n, 6*1) = i°(u) - i^(n) cos 61 . (62) 

In the test-particle limit when the mass of the rocket is negligible (m — )• 0) and A = 0, 
the background flat-space coordinates are obtained by the transformation ()17p . Using 
the relations (jlSp and ()29p this explicitly becomes 



Z° = z°(u) + 



p{u,6i 



= z\u) + —— cos 01, (63) 



1 ^ 
Z = — — — sin 9i cos 82 , 



i-l 

r 



sin ^1 JJ^ sin cos 9i , i = 3, . . . , D — 1 . 



where, as in (f30]l . we define 9d-i = 0. Obviously, at any time the origin r = of the accel- 
erated coordinates of ()6ip coincides with the rocket. 

It now appears to be convenient to introduce the following parametrisation of the com- 
ponents of the velocity u: 

= cosh (/ a{u) dn) , 
z^ = sinh (/ a{u) du) , (64) 

which identically satisfies the condition (j60p . The acceleration vector a has the only nonva- 
nishing components z^ = a{u) z^, z^ = a{u) iP , so that a • w = 0. Moreover, a ■ a = a^{u) 
which means that the function a{u) is exactly the value of instantaneous acceleration of the 
rocket at the time u. Using this parametrisation, the metric function ()62p becomes 

p{u^ 61) = cosh (/ a{u) du) — cos Oi sinh (/ a{u) du) . (65) 

Finally, it is possible to define a new angular coordinate 'd as 

p{u,9i) 
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Straightforward calculation shows that ()66p implies the following interesting identities: 

sin 9i 



sint? 
cos •& 



cosh (/ a{u) du) — cos 9i sinh (/ a{u) du) 
sinh (/ a(n) du) — cos 9i cosh (/ a{u) du) 
cosh (/ Q!(u) du) — cos 01 sinh (/ a{u) du) 



9\ "!? 

cot — = tan — exp (/ a{u) du) , (67) 

(logp),„ = a{u) COST?, 

p^"^ = cosh (/ a(n) du) — cos sinh (/ a{u) du) , 
cos 9i = sinh (/ a{u) du) — cos ■!? cosh (/ a{u) du) , 
p-'^{d9l + sm^9i d6'|) = (di? + a siniJdu)^ + sin^ i?d6'| . 

With these relations, the parametrisation (|63|) of the Minkowski background (when m = = A) 
in terms of the accelerated coordinates becomes 

= 2;'^(u) + r [cosh (/q(u) dn) — cos sinh (/a(u) du) ] , 
Z^ = 2;"'^(u) + r [sinh (/a(n) dti) — cos?? cosh (/q(u) dn) ] , (68) 
= r sin-i? cos 92 , 

i-l 

= r sin sin 9j^ cos ^j, i = 3, . . . , D — 1 . 

i=2 

The complete metric ()6ip takes the form 

D-2 i-l 

ds^ = r2(dT? + a{u) sini?dn)^ + sin^ (d^^ + ^ ( ]J sin^ ej)d02^ 

i=3 j=2 

I , r. 2m(u) 2A 9\ , 9 . s 

-2 d^x dr - (l - 2 r a(n) cos ^ - - r') du', (69) 

which can be rewritten as 

2m(n) 2A 



ds' = -{l 



„ „ — — 2 a(n) r cos ■!? — a^(n) sin^ 7?) dti^ 

(L) — 2)[D — 1) / 

2 du dr + 2 a{u) sin du d'd 



. D-2 i-l ^ 

f di92 + sin^ (d + 5Z ( n ^^f) ) • (70) 

^ i=3 i=2 ^ 

This is a generalisation, to any D > 4, of the standard metric form ([8]) of the Kinnersley 
photon rocket [T][^][71llU| accelerating arbitrarily in a single spatial direction. In the absence of 
acceleration, a{u) = 0, it reduces to the higher-dimensional Vaidya-(anti-)de Sitter solution 
which describes a spherically symmetric fixed source at r = with a varying mass determined 
by m{u). 

Alternatively, the metric for the photon rocket accelerating in a single spatial direction 
(now Od-i) can be written in the Cartesian- like coordinates ()50p . ()5ip with P simplified to 
P = (i° — z^~^) + j{z'^ + z^~^) 5ij x^xK Following (f6ij) . this becomes 

P = exp (— / a{u) du) + \ exp (/ a{u) du) 6ij x'^x^ . (71) 
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6 Explicit examples of motion of a photon rocket 



It will now be illustrative to investigate in more detail some particular examples which 
describe accelerating or decelerating photon rockets along straight or curved trajectories. In 
all the cases we will identify and study the particular situation in which the minimal amount 
of radiation is emitted, i.e., the mass decrease of the rocket (its "energy consumption") for 
a given flight is minimised. 

6.1 Straight flight 

The simplest situation arises when the photon rocket moves in a single spatial direction. It 
does not perform any manoeuvres, except for changing its speed. As described in section [5l 
in such a case the timelike worldline z°'{u) of the rocket is fully described by ()59p . i.e., by 
the two functions z^{u) and z^{u). With the parametrisation of the corresponding velocity 
u given by ([Mj) . 

z^{u) = cosh {J a{u)du) , i^(n) = sinh (/a(u) du) , (72) 

where a{u) is the acceleration of the rocket as a function of the time u, the metric takes the 
form (|70p. It follows from relations ()68p that the angle introduced in ()66p . takes the value 
I? = behind the rocket while -d = ti in front of the rocket (considering a > 0). 

The associated field of emitted photons has the form ()44p where the profile n is given by 
([l5|). Using the relation (logP)^„ = [\ogp)^u = a{u) cos??, see (f37|) and (fBTp . 

D — 2 

n'^{u,'d) = — [-m„(n) + {D - 1) a{u) m{u) cos '&]. (73) 

Svr ' 

Obviously, the maximum of radiation is emitted directly behind the rocket (for = 0) and 
minimum in the direction in front of it (for = vr). In fact, it is possible to rewrite the 
expression ([75]) as 

n'^{u,'d) = ~ am{l + cos'd) - ^^^^ [m ^ + {D - 1) am]. (74) 

Svr Svr 

The first term is always positive, and vanishes for -i? = vr. Therefore, the second term must 
also be nonnegative. An optimised situation occurs when the second term vanishes: in such 
a case the rocket emits no photons in the direction of its motion, it only emits "backwards" . 
This condition yields the explicit relation 

m{u) = mo exp {—{D — 1) / a{u) du). (75) 

The mass of such photon rocket decreases exponentially from its initial value mg, and the 
corresponding radiation pattern (|74p is explicitly given as 

n^iu, 1^) = ^ ^ ( - m^u) cos^ | , (76) 

where 

— m^u{u) = {D — I) ruQ a{u) exp (^—{D — I) J a{u) du^ . (77) 

At a given time u, the angular dependence of the photon field on i} and 6i is plotted as a 
spherical polar diagram in figure El 

In particular, for a photon rocket which moves with a constant acceleration a = const., 
it is possible to calculate its motion and the total decrease of its mass explicitly. Integrating 
(|72p for the initial conditions z^{0) = = z^{0) we obtain 

z^ (u) = — smh.{a u) , z^(ti) = — [cosh(Qn) — 1] . (78) 
a a 
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Figure 2: Section through the radiation pattern of the photon field (1761) emitted by the rocket 
which accelerates along a straight line (moving here to the right), as described by the axisym- 
metric metric (170|) . On the right we plot the same pattern with the trivial angular coordinates 6i 
{i = 2, . . . , D — 2) suppressed. 

In view of the transformation ()63p . these two functions give exactly the time elapsed T = 
= z^{u) and the distance L = = z^{u) of the rocket (located at r = 0) with respect to 
the background inertial frame connected to "the Earth" . Recall that u is related to the proper 
time of the photon rocket. We have thus recovered the well-known relations for (integrated) 
time dilation and length contraction which are valid for a uniformly accelerated objects in 
special relativity. The corresponding total decrease of the mass of the rocket is given by 



Am = mo - m{u) = mo(l - g-^-^"^)""), 

where, considering ([75|) . 

au = arcsinh (a T) = arccosh (1 + a L) . 
The plot of Am/mo as a function of aT for D = 4 is shown in figure [3j 



(79) 



(80) 



Am 



1 








2 



aT 




shows the relative mass decrease of the rocket with the inertial 



time measured "on the Earth" . 



Asymptotically, for large values of u, T and L, the function ()79p reduces to 



Am 



1 - (2a T) 



l-D 



1 - {2a L) 



l-D 



(81) 



which demonstrates that the total loss of mass Am approaches the initial mass mo as 
and L'P, where p = 1 — D < 0. For higher dimensions D, the mass decrease is faster, which 
makes the travel of photon rockets in higher dimensions more demanding. 
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It is also useful to express the mass function ()79p of the rocket in terms of its actual 
speed V with respect to the inertial frame, v = / = tanh{a u) which implies e"" = 
^{l + v)/{l-v): 

mo \l + v) 

This simple formula gives the mass m{v) of the photon rocket after it was uniformly accel- 
erated from rest (v = 0) to the speed v. 

Another explicitly solvable model of a photon rocket moving in a single direction arises if 
the acceleration function takes the form a{u) = tanhu for u >0. In such a case, which may 
describe specific initial phase of acceleration of the rocket, the relation ()72p can be integrated 

to 1111 

fcosh7x+-— — ), i\u) = -(coshu-——), (83) 



2 V cosh M / 2 V cosh u 



so that 



1 vr 
T = z^ (u) = — sinh u + arctan e" , 

1 vr 
L = z^{u) = — sinhti — arctan e" + — . (84) 

Compared to the case of a uniform acceleration a = const., as described by (j78p . it follows 
that for large values of the time u the inertial time T and the distance L traveled are 
asymptotically half of those corresponding to a = 1. The total loss of mass of the rocket is 
now 

Am = mo(l - cosh^^^u), (85) 
which has the same asymptotic behaviour as (jSip for a = 1. 



6.2 Circular trajectory 

Let us also study special non-geodetic flight of the photon rocket, namely a circular motion. 
We assume that the rocket in D dimensions moves along a circle of radius a with a constant 
angular velocity oj, 



z\u) = VT+aho^u, 
z\u) = ... = z^-3(n) =0, 
z^~'^{u) = a cos{uju) , 

z^-^{u) = a sm{uju) . (86) 

Here w is the angular velocity of a photon rocket with respect to the time u. Due to the 
time dilation/length contraction, the constant speed v of the rocket on its circular motion, 
measured with respect to the inertial background frame in the center, is not aoj but 

auj V 

< 1 , i.e., au = 1^ ^ . (87) 



The velocity parameter au may thus take any value, and w — )• 1 as acj — )• oo. 
Explicit exact metric which describes such motion of the photon rocket is 

ds^ = ^ ^ ( n si^^ ^i) '^^i - 2 du dr 



P 1=1 j=i 



2m(u) 2A 2 



l-2r(logp),„--^-^^-^^^^r )dn 
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where 



Z)-3 



p(n, 6j, (p) = \/l + a^a;^ ~ ( IT si'^^i) sin( 



(89) 



ei,92,...,BD-3 G [0,vr] and = 6*0-2 G [0,27r), see ([52]) and (f53|). 

The corresponding radiation pattern of emitted photons is given by Since ( log P)^^ 
(logp)^„ , see (f37|l . this becomes 



in which we introduced 

f{u,9j,4>) 



D 



Stt 



+ {D — 1) auj^m f{u, Oj, (f))], 



n 



f=i'sin^, 



cos((/> — ti) 



D-3 



sin( 



sin(</) — uju) 



(90) 



(91) 



In D = 4 this simphfies to ^ IljLi'^ sin^j^ = sin^i = sin0, and the angular part of the metric 
dSHD reduces to (d6'2 + sin^ ddcj)'^). 

The function / identically vanishes whenever 9j = and 9j = vr which, due to the def- 
inition ([29]) . describe directions perpendicular to the plane of the circular trajectory ([86]) 



of the photon rocket. Also, / = if (j) — cou = ^ or (p — oju = For fixed (j) and u, the 



of motion. If 9 

s\n((h — UJU 



maximum of the function / occurs at 9j = ^ for all j = 1, . . . , — 3, that is in the plane 
J — \ for all j, and u is fixed, the extremes of / are given by the condition 
acj/Vl + cfioj"^ which implies ian{(j) — uu) = zbacj. For (p — uju < ^ the upper 
sign applies and the corresponding direction gives the maximal value of /, namely / = +1. 
Contrary, for (p — uju > ^ the lower sign gives the direction in which there is a minimal value, 
/ = —1, in the radiation pattern. A typical plot of the function f{u = 0,9j,(p) for several 
constant values of 9j is shown in figure HI 

/ 

+ 1- 




ure 4: The function f{u,6j,(p) for u = and several constant values of 6j between and 
we set a = 1 = u. 

Since f{u, 9j, (p) G [—1, +1] and the minimal value / = —1 is attained, to satisfy the condition 



n {u,9j,(j)) > everywhere, we rewrite (j90p as 



n^(ti, 9j, I 



(Z?-l)(Z)-2,_^, , 



■iTT 



auj m[l + f{u, 9j, i 



D 



Svr 



■[m^u + {D -l)auj'^m]. (92) 



Similarly as in (|74p . the first term is nonnegative and vanishes at the minimum of /. An 
optimised flight of the photon rocket, which minimises the mass decrease, thus occurs if the 
second term in (|92p vanishes, i.e., when 



m{u) = mo exp { — {D — 1) auj^ u). 



(93) 



The mass of the photon rocket then decreases exponentially. The coefficient auj"^ is, in fact, 
the classical centrifugal acceleration (notice the analogy with expression (j75p for a straight 
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trajectory, in particular when the acceleration a is constant). The corresponding radiation 
pattern is 

r?{u,Qj,4>) = £^{-m,u)[l + fiu,ej,cP)]. (94) 

As an illustration, assuming D = 4, the angular dependence of the emitted photon field 
on 6 and (for a given value of the time u) is plotted as a spherical polar diagram in figure [5j 




Figure 5: Radiation pattern n^{9,(j)) of the photon field emitted (at u = 0) by the rocket 
moving in four dimensions along a circle of radius a = 1 with the angular velocity u = 1 (left). 
Equatorial section 6* = f through the diagram, where v?{9,(j)) reaches extreme values (right). 




Figure 6: Radiation patterns n^(^, 0) of the photon field emitted by the rockets which move in 
four dimensions along circles of different radii a with the same proper angular velocity uj (here 
oj = 1 and It = 0). For each a the curves plotted indicate the value of n'^{6 = const., (p) in a polar 
graph, with the outer curve corresponding to the equatorial section ^ = f (cf. figure |5] for the case 
a = 1). The asymmetry of the patterns grows for large a and thus velocities v. 

Let us recall that, according to dM]), the radiation density p decreases as p = r'^~^v?{6j, </>), 
where r is the distance from the rocket, namely an affine parameter along null geodesies on 
u = const, generated hy k = dr- 
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In figure [6] we plot the radiation patterns (|94p emitted at u = by the rockets which 
move along circles of different radii a, keeping their proper angular velocity u fixed {uj = 1). 
The curves shown correspond to 9 = const., i.e., for a given a and 6 the value of n'^{9, (j)) is 
plotted radially with (p being a standard polar angle. For larger o — and thus larger circular 
velocity v given by ()87p — the radiation patters become more distorted. 

Due to the presence of the argument (j) — lou in ()9ip and subsequent expressions, it is 
clear that the radiation pattern rotates along the circular trajectory ()86p with the proper 
angular velocity u. This is shown in figure [7] for the choice auj = 1 and L) = 4. Moreover, the 
radiation of photons decreases exponentially with the factor exp(— 3 a(^^ n), see ()94p . ()93p . In 
the proper scaling, the smallest pattern in figure [7| (on the left) should, in fact, be e^'^-times 
smaller than the largest one (on the right). Physically, since the mass of the rocket decreases 
exponentially, an exponentially decreasing "reactive force" from the emission of photons is 
sufficient to keep it on a circular motion with constant speed v. 



Figure 7: Radiation pattern n^{6,(f)) of the photon field emitted by the rocket at several places 
during its circular motion around the centre {auj = 1, D = 4). The dashed lines indicate directions 
in which the radiation is maximal and minimal. 

We can also evaluate the mass function ([93]) after the "U-turn" depicted in figure [71 using 
([87]) and the final time uu = n, as 



Let us now compare this expression, which gives the mass loss after the circular manoeuvre 
during which the speed v of the photon rocket is kept constant while the direction of its 
flight is reverted, with the analogous relation corresponding to a constant deceleration q, 
followed by the same acceleration, along a straight line. Specifically, we employ the relation 
()82p which gives the mass m{v) of the rocket after its uniform acceleration along a straight 
line from v = to the speed v. However, this has to be combined with the initial straight 
flight during which the photon rocket decelerates from the velocity —v to zero. By extending 
relations ([78p . ([79]) and ([82]) to negative values of u and v we obtain 





(95) 




(96) 
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Since 



TTV 



< log 



1 



for any v £ (0, 1), the mass function ()96p is greater than 
functions in figure [HI In particular, for a smah speed v, 



(97) 

see also the plot of these 



m{v) 
m{—v) 



1-2{D -l)v , while 



miy) 
mo 



1 -vrp -l)v. 



(98) 



We can thus conclude that in any dimension D the deceleration from the velocity v to zero 
and subsequent "backward" acceleration to the same velocity v in exactly opposite direction 
along a straight line is more efficient than the circular U-turn of the rocket because the total 
mass of the photons emitted (the "energy consumption") is smaller. 



1.0 



0.8 



0.6 



0.4- 



0.2 




Figure 8: Plot of the mass functions (l96i) and ( 195|) . assuming D = A. The former is greater than 
the latter for any speed v G (0, 1). 



7 On motion of test particles and absence of "grav- 
itational aberration" 

In this final section we will briefly comment on geodesies which describe motion of free test 
particles in the exact spacetimes representing the gravitational field of an arbitrarily moving 
photon rocket. 

To investigate such geodesies we have to evaluate the Christoffel symbols V^^ for the 
metric ()42p . With respect to the Minkowski coordinates Z", the metric reads 

9c.p = r]ap + 2VKkp, so that g''^ = - 2V k^^k^ , (99) 

where V{r^u) is given by ()43p . Moreover, using (jl7p and ()12p - ()14p . for the functions m(Z") 
and r{Z°') we obtain the relations 

= (1 + rfc^r ) A;^ - , (100) 

ij'ka) ,^ — "Uaii ~l~ k^ , 

where Za = t/q,^ . Straightforward but somewhat lengthy calculation then yields 
2V 

= — ^" VfS-y + Vr Z'^kpk^ + {2VVr " Vu) k'^kpk^ 

+ (— - Vr) {k'^zpk^ + k'^i^kfi - (1 + rk^r) k'^kpk^) . (101) 
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Considering = —k^ = p ^ , k^ = ki = p ^n* where p is given by ()22p and n = n*(9j is the 
spatial unit vector (see section ED , all components of F^^ can be written explicitly. 

It is also convenient to introduce a spatial velocity vector v = v^di of the rocket with 
respect to the Minkowski background frame by f * = i* (u) / (u) . Consequently, 

i°(u)=7, i*(u)=7v% where 7 = , (102) 

^/l — V ■ V 

is the standard Lorentz factor corresponding to the velocity v of the rocket at V (i.e., at its 
proper time u), see figure [H The function p then simplifies to 

p = j{l-n-v). (103) 

Let us recall that n is the unit vector connecting the (spatial position of) events V and Q. 
In fact, the spatial distance R between V = {z^{u), . . . , z^~^{u)} and Q = {Z^, . . . , Z^^^} 
isR= ^(Zi - zi(ti))2 + • • • + (Z^-i - z^-'^iu))^, which using ([H]), (HH) becomes R = rk°. 
It follows that 

r = pR. (104) 

Prom (fTUT]) . (HHD we thus obtain 

_ / 2m 2Ai? \/ (I)-3)m 2Ai?2 

J- 00 pD-2„D-i m _ 9Vn _ 1 w ^ 




RD-2pD-i p_2)(D-l)M RD-^pD-i {D-2){D-1) 



{3-D)m 2AR 



RD-2pD^l (D-2){D-1 

/(£-l)(27p-l)m (Z)-l)mfc^z'^ + m ^ , 

+ ^ RD-2pD+l J^D^ ) " ■ 

Due to the geodesic equation, the coefficient —Tqq gives the "spatial acceleration" (in New- 
tonian terminology) of a test particle which is at rest at a given point Q. Obviously, for 
large distances R of the test particle from the photon rocket, the dominant contribution 
to its acceleration arises from the cosmological constant A which represents global isotropic 
expansion of the spacetime when A > 0. On the other hand, for small R the contribution of 
the A-terms is negligible. 

Setting A = in (jlOSp and employing relations (jl02p . (jlOSp we obtain 



-L n 



00 



RD-2^D~l {l_n- V)D+^ 



[{D - 2,) - 2{D - ?,)n ■ V - 2{n ■ vf 



+{D -l)v.v]v}-{D-2.){l-n-v)v'- , ..ad^s " 



(D — l)m k„z" + m ■ 

RD-3jD(^l_n.v] 



D 



(106) 



The last term, proportional to the acceleration z'^ and mass decrease m of the rocket, rep- 
resents the radiative part of the gravitational field which behaves as oc R^-^. Close to the 
photon rocket this can be neglected with respect to the first "Newtonian force"-term oc R?~^ . 
Notice that such Newtonian-like gravitational acceleration is oriented along 

{D - 3)[(n^ - v') - {2n' - v') n ■ v + . . .] , (107) 

which is the spatial direction toward the "instantaneous" position of the rocket, extrapolated 
from its "retarded" position given simply by n*. This demonstrates the absence of "gravita- 
tional aberration" in such systems which are explicit exact solutions of Einstein's equations 
in any dimension D. For a thorough discussion of this issue in L> = 4 see the Carlip work 
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8 Conclusions 



We analysed in detail some properties of the class of exact spacetimes which represent Kin- 
nersley's photon rockets moving arbitrarily in any dimension D. These solutions contain 
pure-radiation photon field and admit a cosmological constant A. They are of algebraic 
type D and can be written either in the Kerr-Schild form ()42p . ()43p or more explicitly in the 
Robinson-Trautman forms ([50]) . ([5T]) and ([52]) . ([53]) . 

In section [3] we first systematically reviewed general Newman-Unti [ (|2ip . (|22p ]. spherical 
[ (|2D . ([25]), or dSZl), ([28]), or (IMl), dM])] and Cartesian-like [([MD-dMl)] background coordinates 
suitable for the description of an arbitrarily moving test source in flat Minkowski space of 
dimension D. 

The complete Kinnersley solution in various metric forms was presented and discussed 
in subsequent section HI A special case of photon rockets accelerating in a single spatial 
direction is contained in section [SJ 

We discussed important particular trajectories of the rockets, namely a straight flight (in 
subsection 16. ip and a circular motion (in subsection I6.2p . including the corresponding radi- 
ation patterns of the photon field and the mass loss formulae. For example, we showed that 
the straight deceleration followed by the backward acceleration is a more efficient manoeuvre 
then the circular U-turn since the total mass of the emitted photons (and thus the energy 
consumption) is smaller in the former case. 

In the final section [7] we derived the Christoffel symbols for general Kinnersley spacetimes 
which are crucial for discussion of geodesies. In particular, we demonstrated that, in any 
dimension D, there is no "gravitational aberration" effect. 

To conclude, let us remark that in our contribution (as well as in previous works [2lll0j) 
the Kinnersley solution has been used as an exact model of an accelerating photon rocket. 
However, in D = 4, it may also be considered as a simple relativistic model of non-geodetically 
moving astronomical object in the Solar System (such as a small particle or an asteroid) 
which accelerates due to a specific anisotropic thermal absorption/emission of photons, i.e., 
to simulate the Poynting-Robertson effect |21] or the Yarkovsky effect |22j . 
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